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ABSTRACT 

In the Excursion Set approach, halo abundances and clustering are closely related. 
This relation is exploited in many modern methods which seek to constrain cosmolog- 
ical parameters on the basis of the observed spatial distribution of clusters. However, 
to obtain analytic expressions for these quantities, most Excursion Set based predic- 
tions ignore the fact that, although different /c-modes in the initial Gaussian field are 
uncorrelated, this is not true in real space: the values of the density field at a given 
spatial position, when smoothed on different real-space scales, are correlated in a non- 
trivial way. We show that when the excursion set approach is extended to include such 
correlations, then one must be careful to account for the fact that the associated pre- 
diction for halo bias is explicitly a real-space quantity. Therefore, care must be taken 
when comparing the predictions of this approach with measurements in simulations, 
which are typically made in Fourier-space. We show how to correct for this effect, and 
demonstrate that ignorance of this effect in recent analyses of halo bias has led to 
incorrect conclusions and biased constraints. 
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1 INTRODUCTION 

The abundance of massive galaxy clusters is expected to 
be related to the abundance of sufficiently overdense iso- 
lated regions in the initial conditions (Press & Schechter 
1974). In the excursion set approach (Epstein 1983; Bond et 
al. 1991; Lacey & Cole 1993), the abundance dn(m)/dm of 
such overdense 'clouds' containing mass m that are not con- 
tained within larger overdense clouds is estimated by map- 
ping the problem to one involving random walks. If f{s)ds 
denotes the first crossing distribution of a barrier of a cer- 
tain prescribed height by random walks, i.e., the fraction of 
walks which have crossed a barrier B{s) in the 'time' inter- 
val (s, s + ds), without having crossed it before, then the 
Excursion Set ansatz is that 



m dn{m) 
p dm 



dm = /(s) ds. 



(1) 



where p is the background density. In this approach, the bar- 
rier height and its dependence on s is set by the competition 
between gravitational collapse and the cosmological expan- 
sion history (Sheth, Mo & Tormen 2001), and the relation 
between random walk time s and halo mass m is set by the 
statistics of the initial fluctuation field. 
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To obtain analytic expressions for the first crossing dis- 
tribution, Bond et al. (1991), and most analyses which fol- 
lowed, assume that the density field smoothed on one scale 
is trivially correlated with that on another scale. Namely, 
if one plots the overdensity (5 as a function of smoothing 
scale s, then this resembles a random walk in 'time' s - the 
usual assumption is that successive steps in the walk are 
independent of the previous ones. This assumption is only 
correct for a Gaussian field smoothed with a series of fil- 
ters that are sharp in fc-space; it is incorrect for any other, 
more physically reasonable, set of smoothing filters such as 
a Gaussian or a TopHat in real-space (Bond et al. 1991). 
For walks with uncorrelated steps, the first crossing problem 
can be solved exactly for barriers whose height is a constant 
(Bond et al. 1991) or a linear (Sheth 1998) function of ran- 
dom walk time. In particular, the first crossing distribution 
of -B(s) = (5c (1 + Ps/5l), for some slope P, is 

5(0) 



Sfu{s) = 



-S^(s)/2s 



(2) 



where the subscript u is to remind us that this is for walks 
with uncorrelated steps. Reasonably accurate analytic ap- 
proximations to the exact solution for more general barri- 
ers are also available (Sheth & Tormen 2002; Lam & Sheth 
2009). 

In this approach, the question of halo bias reduces to the 
problem of estimating the first crossing distribution when 
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the walk starts from some prescribed height on some pre- 
scribed scale {S, S), rather than from the origin (Mo & White 
1996; Sheth & Tormen 1999). The halo bias factors bk are 
defined by expanding the conditional distribution in a Tay- 
lor series around 5 — and S — and ignoring all terms of 
order S/s: 



f{s\S,S) 
fis) 



(3) 



The coefficient of the first term in this series, the linear bias 
factor bi, is the one of most interest. 

For a barrier of constant height B{s) — Sc and walks 
with uncorrelated steps, except for a shift of origin, the con- 
ditional walk follows the same statistics as the unconditioned 
one. Namely, one simply sets Sc ^ 5c — 5 and s — >■ s — S in 
the expression for the first crossing distribution. Therefore, 
in the S ^ s and 5 ^ 5c limits, instead of expanding the 
conditional distribution in a Taylor series (as in equation[3]), 
the halo bias factors may equally well have been obtained by 
differentiating the unconditional distribution with respect to 
5c. This remains true for more general barriers: if we use 5c to 
denote the barrier height on very large scales, i.e., 5c — B{0), 
then 



bu = 



d In /„ (s) 
d5c 



+ 13-1 



(4) 



where the derivative with respect to 5c actually means a 
derivative with respect to S(0), and the final expression on 
the right assumes B{s) = 5c (1 + l3s/5^). (When /3 = 0, then 
our bu reduces correctly to the bias formula given by Cole & 
Kaiser 1989 and Mo & White 1996 for a constant barrier.) 
It is in this sense that a model for halo abundances carries 
with it a model for halo bias. And it is this close relation 
between halo abundance and bias which vastly simplifies 
analyses which seek to constrain cosmological parameters 
from cluster abundances and clustering. 

Recent work has shown that this close relationship is 
not accurate at ten percent precision (Manera et al. 2010). 
This is one of the reasons why there is renewed interest in go- 
ing beyond the assumption of uncorrelated steps. Peacock 
& Heavens (1990) provided an analytic approximation to 
the first crossing distribution of a barrier of constant height, 
when the steps are correlated because of a Gaussian smooth- 
ing filter. Recently, Maggiore & Riotto (2010) have provided 
another approximate solution to this problem, which is bet- 
ter motivated for a real-space TopHat. As Paranjape et al. 
(2011) discuss, the two approaches may be thought of as ex- 
panding around the two opposite limits of completely cor- 
related and completely stochastic walks, respectively. Since 
Monte Carlo solutions of the exact answer show that the first 
crossing distribution for a TopHat is always within about 
ten percent that for a Gaussian (e.g. Bond et al. 1991), the 
question arises as to which approach is the more accurate. 
Paranjape et al. showed that it is the older Peacock-Heavens 
approximation which is the more accurate. 

Recently, Ma et al. (2011) have extended the Maggiore- 
Riotto approach to derive expressions for halo bias. They 
find that this extension appears to predict larger bias fac- 
tors than for the case of uncorrelated steps (i.e., smoothing 
filters that are sharp in fe-space). Since the Peacock-Heavens 
approximation to the first crossing distribution is the more 
accurate, the main goal of the present paper is to compute 



the associated prediction for halo bias, to see if it too leads 
to enhanced bias factors. We do this in Section (2] finding 
that it does. 

However, our analysis highlights an important property 
of smoothing with a sharp-A: filter. Namely, for this filter, 
halo bias does not depend on whether one measures it in 
real-space, by measuring the ratio of halo counts-in-cells 
conditioned to have a certain overdensity to the uncondi- 
tioned average, or from the ratio of the halo-halo or halo- 
mass power spectrum to that of the mass. The excursion 
set prediction is explicitly for the former, whereas measure- 
ments in simulations are almost always of the latter. We 
show that when one uses anything other than a sharp k- 
space filter then the bias measured from halo counts-in-cells 
differs from that measured using power spectra. This is the 
physical origin of the enhanced bias associated with TopHat 
or Gaussian smoothing filters. A final section discusses some 
implications. 



2 HALO BIAS AND CONSTRAINED WALKS 

For what follows, it will be useful to define 

^m^j^^'^k^'W\kRl (5) 

where Ps{k) is the power spectrum of initial density fluc- 
tuations (linearly extrapolated to present epoch), and W is 
a smoothing fllter. The quantity o-qIR) measures the vari- 
ance in the field on scale R. We will often denote this as 
s = (Toi 8'iid will only show the _R-dependence when omit- 
ting it would have led to ambiguity. In addition, we will 
also make frequent use of the correlation scale defined by 
Peacock & Heavens (1990): 



27rln(2) ^77(1-7'), 



where 



(T0O"2 



(6) 



(7) 



For Gaussian smoothing filters and Ps{k) oc fc", which we 
will use to illustrate our results, (j| oc jj-3-"-2i g^jj^j ^2 _ 

{3 + n)/{5 + n), making Fph = 2 7rln(2) ^(n-f 3)/2. 



2.1 Constrained walks with correlated steps 

The analysis below is restricted to barriers B{s) for which 
B{s)/y/s decreases monotonically. Define 



P[s\5,S) = - 



where 



Via = 



B{s) 



and 



dk k^Ps{k) 



k 27r2 



W{kRs)W{kRs) 



(8) 



(9) 



(10) 



Let f{s\5, S) be the first crossing distribution of the barrier 
B{s) for walks with correlated steps that are constrained to 
pass through some (5, S) with S < s, while also remaining 
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below the barrier B{S') on scales S" < 5*. Then a good 
approximation to f{s\S,S) is, 



f{s\S,S) 



d 



P{s\5, S)E{s\5, S) 



where 



E{s\5, S) = exp 



As' 



\aP{s'\5,S) 



(11) 



(12) 



and FpH was defined in equation ((Bj (Paranjape et al. 2011). 

The unconditional first crossing distribution f{s) is 
given by setting r = in the expression for uiq: 



P(s) = i [l + erf (i3(s)/^)] , 
E{s)=exp( [ 3^ InP(s') 



sTf 



sfis) 



d 



-s—{P{s)E{s)) 



where 



s/c(s) = - 



E{s) 



dP 



sUs)-—P{s)\nP{s) 

1 PH 



(13) 
(14) 

(15) 



ain[B(s)/yi] B{s) 



91n s 



91ns 



B(0) 



sfu{s). (16) 



The limit Ffh — >■ oo is what Paranjape et al. (2011) called 
completely correlated walks. In this limit E ^ 1: the quan- 
tity /c(s) is therefore the first crossing distribution for com- 
pletely correlated walks. Paranjape et al. showed that this 
limit provides an excellent description of the unconditional 
first crossing distribution of barriers which decrease steeply 
with s. In particular, for such barriers, this limit is a better 
description of /(s) than if one uses the value for Ffh that 
is given by equation ((6|. This will be important below. 
Note that /c 7^ /u. For barriers of the form 



B{s)^5,{l+I3{s/Slr) 
where we assume a > and /? < 0, 



(17) 



Sfc{s) = 



B{0)' 

^(l-/3(2a-l)(V4^)") 



(18) 



where we have explicitly shown the dependence on B{0) in 
the first line. For a barrier of constant height (/3 — 0), we 
have fc ~ fu/^. For a linear barrier {a = 1), /c = (1 — 
/„/2. 

2.2 Bias as the large scale limit 

Notice that the integral in our expression for r is similar to 
that which defines s = aQ{Rs) and 5* = a^lRs), the only 
difference being that here the two smoothing filters have 
different scales. If we use Sx to denote the value of this 
integral, then r^Ja/S — Sx/S. This quantity depends on 
the form of the smoothing filter. For a sharp-A; filter (the 
one associated with uncorrelated steps) Sx ~ S, so Vio — 
[B{s) — So]/\^ s — S. Thus, uio is related to the unconstrained 
quantity (r = in equation[9ll by a shift of origin. Therefore, 
the first term bi in the expansion of equation Q will be the 
same as 6„ of equation (|4]). 

For Gaussian filtering of a power-law spectrum 



r = [2RsRs/{Rl + i??)]("+''/' and Sx/S = [2/(1 + 
R^/Rg)]^"^^^^^ . In this case, if one wishes to think of uio as 
an effective shift of origin, then this shift is scale dependent. 
However, ff Rs <S:Rs, then Sx/S 2^"+^'>^^: in this limit, 
the Gaussian case is similar to the fe-space one, except that 
the shift of origin is 4 - 2*"+^^''^(5. This shift comes from the 
fact that steps are correlated. As a result, following equa- 
tion (|3]) above will amount to the same as differentiating the 
first crossing distribution with respect to 5c (or more gener- 
ally, with respect to B{0)), only taking care to account for 
the factor of Sx/S 2'^"+^^^^. 

The discussion above applies for other smoothing fil- 
ters too, except that the numerical coefficient 2^"'^^^^'^ will 
change. We will argue shortly that the actual value mat- 
ters little for our purposes, but to complete the discussion, 
we provide some explicit examples of this change. For a fil- 
ter that is a tophat in real space, Sx/S — >■ 5/4,4/3 and 
1 for n = —2,-1 and 0. Thus, in general, the numerical 
coefficient is smaller than when the filter is Gaussian. But 
the important point is that, for walks with correlated steps, 
equation should be replaced by 



bi 



Sx_\ dlnfjs) 
S J dSc ' 



(19) 



where Sx/S ~ 1 depends on the filter and the power spec- 
trum, and /(s) is the unconditional first crossing distribu- 
tion. Note that, to obtain this result, we use equation (|8]) in 
the limit S <ti s, so that we have uu) — >■ (Sc — SSx / S) / ^ 
where Sx/S is constant. Plugging this into equation (|lip . 
expanding to linear order in S, and using equation (|3]) finally 
results in equation ([19}. 



2.3 Bias in real vs Fourier space 

To better appreciate the origin of the extra Sx /S factor, it 
is useful to consider the case of peaks in the initial Gaussian 
fluctuation field. Since the property of being a peak is scale- 
dependent, we must flrst specify the smoothing-scale and 
filter with which the peak was identified. We will use _Rp]j to 
denote this scale. 

Bardeen et al. (1986; hereafter BBKS) provide expres- 
sions for the mean density profile around a peak (their equa- 
tion 7.8 and Appendix D) - as they note, this is essentially 
the same as the cross-correlation function between peaks 
and the dark matter field - and for how the large scale den- 
sity modulates the abundance of peaks (their Appendix E). 
Each of these is explicitly a real-space statement, and so 
BBKS provide different expressions for each case. Note that 
the former is one way in which peak-bias is defined, and the 
latter is close in spirit to the conditional crossing calculation 
in the previous section, where we were interested in how the 
first crossing distribution was modulated by knowing the 
value of the density field on some large smoothing scale. 

However, Desjacques et al. (2010) show that all these 
BBKS expressions are in fact a consequence of the fact that 
peaks may be thought of as being linearly biased tracers of 
the dark matter, with a fe-dependent bias factor: 



5pk(k) = fepk(k) S^{k)WikR 



pk) 



(20) 



The property of being a peak is scale-dependent, which is 
why the smoothing scale appears explicitly in the expres- 
sion for the peak field. And note that the value of the bias 
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factor actually depends on properties of the peak (height 
and curvature), so it depends implicitly on i?pk. That is 
to say, peak-bias is simple in Fourier space, and all the key 
BBKS expressions for real-space quantities (e.g. in their Ap- 
pendices D and E) are simply consequences of multiplying 
the above bias relation with different filter functions before 
Fourier transforming to obtain the real-space expressions. 

In particular, the cross-correlation between peaks (de- 
fined on scale -Rpk) and the mass smoothed on scale Rs is 



For the barrier of equation (|f 7|) . we have 



SpkSs ) = 



27r2 



WikRpy,)W{kRs) (21) 



where we have defined 3$ = Sm{Rs)- In the high-peak limit 
(and in the fc ^ 1 limit), b{k) becomes independent of k. In 
this case 



and so 



<5pk|5s 



SpkSs 



bS^ 



Ss 



S 



Ss. 



(22) 



(23) 



This shows that when the Fourier space bias factor is b, then 
in real space, the bias picks up a factor of Sx/S. 

It is usual to assume that high peaks are not a bad 
model for massive halos. So it is not unreasonable to assume 
that the expression above, with the replacement Spi^ — > 5h is 
appropriate for halo bias. Since the excursion set approach 
identifies 5h with the right hand side of equation (|3]), com- 
parison of equation (|23|) with equation (|19p suggests that 
din f{s) /dSc returns the more fundamental Fourier space 
bias factor, even though /(s) was estimated from a real- 
space calculation; the factor Sx/S arises simply because 
the excursion set calculation returns the cross-correlation 
between halos (regions in the initial conditions which are 
above Sc on scale Rs) and the initial mass fluctuation field 
(smoothed on Rs ^ Rs)- 

In the Appendix, we discuss the consequences of equa- 
tion (|20|l with 5pk — > 5h and _Rpk — >■ Rh, for real-space 
2-point correlation functions (in which the matter field is 
not smoothed on some large scale -Rs). In particular, the 
bias determined from such real-space correlation functions 
will depend on halo mass, with the differences between the 
different measurements being in qualitative agreement with 
results of A''-body simulations (Manera et al. 2010). 



2.4 Halo bias for constant or decreasing barriers 

Having motivated why d In f{s) /dSc is fundamental even for 
walks with correlated steps, we now evaluate it for a few 
special cases. (Recall that the derivative here is actually with 
respect to B{0)). The general expression is 



ain/(s) _ d\nE{s) 



85. 



+ 



din Ms) 



1 dP{l + \nP) 



FPH dSc 



'Ms 



1 - 



1 P InP 



(24) 



where P{s) and /c 
H16p . respectively. 



FpH sMs) . 
were defined in equations (|13p and 



din E{s) 
8P{s) 



1 



8 In Ms) 
85. 



5cFpH Jq 
1 

1 f Sfn{s 



ds 



, Ms') 

Pis') ' 



-sMs) : 



5. \2sfc(s)^ 
X (^l-5,Vs + 2/3(Q-l)(s/5e')""' 

+ l3\2a-l){s/5!)^ 



(25) 



For the constant barrier (/3 = 0), din E{s)/dln5c = 
—2 lnP(s)/rpH. In general, equation (I24|> shows that, as 
FpH oo, 8ln f{s)/85c — > 8ln fc{s)/85c, as it should (the 
term involving the derivative of E is always proportional to 
I/Fph). For linear barriers (a = 1), this limit provides a 
better description of the first crossing distribution for suf- 
ficiently negative /3 (Paranjape et al. 2011). Therefore, we 
might expect this to be true for the bias factor also. In this 



91n/c(s)_ 1 fl-5^/s + l3'^s/5l 



Ps/51 



(26) 



2.5 Comparison with Monte-Carlo solution 

To test our analysis. Figure [1] compares the bias factor we 
measure by direct Monte-Carlo simulation of bi via equa- 
tion ((3}. I.e., we construct the first crossing distribution of 
B{s) = 5c by walks conditioned to pass through some S <^ 5. 
on scale S ^ 5^. Then we divide this by the (Monte-Carlo 
estimate of the) unconditional first crossing distribution, 
and subtract 1. Our Monte-Carlos assumed a power spec- 
trum of the form Ps{k) oc k~^'^ and Gaussian smoothing fil- 
ters, for which Sx/S 2°'^ and FpH = 4.3. The solid curve 
shows the result of inserting equation (|24|) with Fph = 4.3 
in equation p9p . It provides an excellent description of our 
measurements. The dotted curve shows the bias relation for 
uncorrelated walks (i.e. a sharp-fe filter), equation Q, is not- 
icably different, lying well below the measurements at high 
Sc/a. 

The dashed curve shows the corresponding expression 
for the bias from Ma et al. (2011): the right hand side of 
their equation (24) minus 1, times 5.. Their expression has 
one free parameter, k; we set it to 0.35, the value they say 
is appropriate for Gaussian smoothing filters. Notice that, 
although it is not as good a description of our measurements 
as is our solid curve, it does a reasonable job, indicating that 
the technical machinery which went into generating this pre- 
diction has produced a reasonably accurate result. However, 
this same technical machinery has obscured the reason why 
their formula predicts a bias factor that is approximately 
twice as large as that given by equation when 5c/(y 3> 1. 
Ma et al. attribute this difference to what they call stochas- 
ticity in the barrier height - but this cannot be correct, 
because stochasticity is not present in our constant barrier 
Monte-Carlos. 

To make this point more directly, note that on the large 
scales which are appropriate for the bias calculation, our 
analysis shows that the conditional distribution /(s|<5, 5* — )■ 



Bias with correlated steps 5 




-0.5 



0.5 



Figure 1. Large-scale linear bias factor associated with a barrier of constant height 5c. Symbols with error bars show the bias factor 
determined via Monte-Carlo solution of equation ([Sj, for walks with correlated steps because of Gaussian smoothing and Ps{k) Qc k~^ -^. 
I.e., Scbi was determined by considering walks conditioned to pass through 5 <^ (5c on scale S <^ S^, for a range of choices of S. Specifically, 
we show [f(s\S, S)/f{s) - l]/(5/5c) for S = (0.155c)^ and S/Sc bin centers ranging from [-1.8, -0.2]\/S, and from [0.2, 1.8]VS, in steps of 
0.4\/5, each of width 0.2\/S. (For clarity, at each v, symbols for different 5 have been offset slightly.) Dotted curve shows 5cbu = — 1, 
which would be appropriate for walks with uncorrelated steps (equation |4j . Solid curve shows the result of inserting equation II 241 1 with 
FpH = 4.3 in equation (I19I I. Dashed curve shows ^ct'l determined from equation (24) of Ma et al. (2011) with their k = 0.35. 



0) should be very well approximated by the unconditional 
distribution /(s) with a shift of origin by {Sx/S)5. Upon 
noting that Ma et al.'s k is our Sx / S — 1, it is easy to 
show that their equation (A32) for /(s|<5, 5* — >■ 0), is indeed 
what our analysis suggests, except that they have used their 
own expression for f{s), and the result has been expanded to 
lowest order in k. However, this expansion in k has obscured 
the simple, intuitive relation between the conditional and 
unconditional distributions that our analysis exploits (e.g. 
Ma et al. themselves failed to notice it). This also explains 
why their expression for the bias factor does not provide as 
good a description of our Monte-Carlos as does ours; their 
expression assumes n <^ 1 when it is not (this matters more 
than the fact that their expression for /(s) is not as accurate 
as ours). 

The results of the previous section suggest that the main 
difference between the solid and dotted curves (for Gaussian 
and sharp-fc filtered walks) is almost entirely due to the fact 
that the excursion set measurement yields a real-space bias 
factor which carries an additional factor of Sx / S: for our 
Gaussian smoothed walks, this factor is 2"'^, whereas it is 
unity for the sharp-fe filter associated with equation (|4]). To 
show this explicitly, Figure[2]shows the result of dividing our 
measurement by Sx / S so as to compare it with the deriva- 



tive of the unconstrained first crossing distribution, which 
we compute analytically using equation H24|) . and show as a 
solid curve. Notice that now the measurements (and our 
solid curve) are much closer to the dotted curve (which 
shows equation |4)| . The differences between the solid and 
dotted curves in this plot are entirely a consequence of the 
fact that the underlying first crossing distributions are dif- 
ferent; they are a much fairer depiction of the effect of cor- 
related steps on the large scale bias. The graphic differences 
between Figures [T] and [5] demonstrate explicitly that cor- 
rectly accounting for the factor of Sx / S is essential; failure 
to do so leads to incorrect conclusions about the nature of 
the relationship between the derivative of the first crossing 
distribution and the large scale bias. 

Figure [3] shows a similar analysis of the bias associ- 
ated with a linearly decreasing barrier (we set a = 1 and 
P = —1). In this case, hu of equation Q lies well below 
the Monte-Carlo solution, whereas b\ of equation (|24|l lies 
above it, at least at small 5c/o". However, the FpH — >■ oo limit 
of equation (|24|l . namely din fc{s) /dSc from equation (|26|) . 
provides an excellent description. This is very reassuring be- 
cause, as we remarked previously, /c(s) itself (the Fph — >■ oo 
limit of the first crossing distribution) also provides a better 
description of the Monte-Carloed first crossing distribution. 
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Figure 2. Same as previous Figure, except that now we have 
rescaled all quantities by a factor of S/ Sx ■ For the measurements 
and the solid curve, this factor is l/2'''^; for the dotted curve it 
is 1. 

We note again that accounting for the factor of Sx /S was 
cruciaL 



3 DISCUSSION 

We have presented an analysis of the large scale bias factor 
associated with the Excursion Set approach for smoothing 
filters which give rise to walks with correlated steps. This 
required careful consideration of the differences between real 
and Fourier space bias - something that has been all but 
ignored to date since, for a sharp-fe filter, these two are the 
same. For other filters, the real-space bias will generically 
be different. 

In particular, our analysis shows that the first crossing 
distribution f{s\S,S), of walks that are constrained to pass 
through 5 on scale S should, as 5* —> 0), be very well approx- 
imated by the unconditional distribution f{s) with a shift 
of origin by {Sx / S) 5 (equations [8] and |9} . Thus, memory 
effects in the walks manifest in two conceptually different 
ways: the fact that Sx/S 1, and that the functional form 
of /(s) itself is altered. We argued that the former is sim- 
ply associated with smoothing effects associated with the 
precise way in which the bias factor was defined; the latter 
is the more fundamental difference. The usual procedure of 
differentiating the mass function with respect to the height 
of the walk on large scales provides a good description of 
the latter, which we identified with the Fourier-space bias 
factor; it underestimates the real-space bias estimated by 
the cross correlation between halos and the mass (smoothed 
on large scales), by a factor Sx/S that can be as large as 2 
(equation (|23[) : this factor is slightly smaller for ACDM-like 
power spectra), while slightly overestimating the bias from 
real-space correlation functions (equation (|X3))). 

Comparison with measurements in Monte Carlo simu- 
lations of constrained random walks showed good quanti- 
tative agreement with our excursion set predictions of the 




lg(i/^ = 6l/a^) 

Figure 3. Same as Figure [2] except now the barrier is B{s) = 
(5c(l + /3s/<52). Symbols with error bars show the bias factor 
determined via Monte-Carlo solution of equation JSjl; dotted 
curve shows the prediction for walks with uncorrelated steps 
(from equation |4]l ; solid curve shows equation I I24I I , and dashed 
curve shows din fc{s)/dlnSc (from equation l|26| l'). which is the 
Fpjj — > oo limit of equation II 241 1. 



real-space bias. In particular, the effect of the Sx /S term 
(the first of our two effects) can be large (compare the solid 
and dotted curves in Figure [T|. The fact that the first cross- 
ing distribution / (equation (|15|) ) is not the same as the one 
for uncorrelated steps fu (equation ([2])) also matters, though 
it is subdominant (solid and dotted curves in Figure [2] are 
different, but this difference is smaller than in Figure [T|. 
However, this second effect is the only one which matters 
for the Fourier space bias. The agreement between the solid 
curves and the measurements in Figures [T] and [2] suggest 
that any additional effects missed by our analysis must be 
small. 

Our results highlight the fact that, when bias is sim- 
ple in Fourier space, then, once one goes beyond sharp-fc 
smoothing, almost all measures of the large-scale bias factor 
will differ from one another, but these differences can be un- 
derstood simply as a result of applying different smoothing 
filters to the same underlying Fourier space bias. Differences 
between real and Fourier space measures of the large scale 
bias factor are now being seen in simulations (e.g. Manera 
et al. 2010), so future comparison of these bias factors with 
simulations should be careful to specify which measure of 
the bias is actually being studied. 

We noted that our excursion set calculation of the dif- 
ference between real and Fourier space bias is consistent 
with the differences between real and Fourier space bias for 
peaks reported by Desjacques et al. (2010). The peaks cal- 
culation is remarkable because there, the bias factors are 
fc-dependent. We believe that one of the virtues of our excur- 
sion set analysis is that it shows clearly how the excursion set 
approach too can lead to fc-dependent bias. Namely, if halos 
of mass m formed from peaks in the initial density field, then 
one would replace equation ((9)l with the expression which 
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describes the density run around peaks rather than random 
positions of height (5c/cr(m). This would yield real-space bias 
factors which, following Desjacques et al. (2010), could then 
be interpreted in terms of fc-dependent bias in Fourier space, 
in effect providing the first calculation of fc-dependent bias 
from excursion set theory. 
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W{kRh) = e 

Cmm(r; Rc) 



we have 
5m(0)(5,„(r)) 
dfc k'-^Ps{k) 



(,hm.{r;Rc,Rh) 



k 27r2 

ar/Rc) , 

dfc fc='Pi(fe) 
k 



jo(fcr) 



(Al) 



= b 



Rh 



2tt2 

2 \ -2 



jo{kr 



1 + -^) ^r/^Rl+R^ 



(A2) 



with the same function ^{x) appearing in both expressions, 
but with different arguments. For Ru ^ -Rc we can then 
write 



1- 2 + 



1 ding 

2 dins 



Rc 



+ . 



(A3) 



where x = r/Rc- For separations r such that g(a;) oc 
this bias from correlation functions therefore picks 
up a relative decrement of |Afe|/6 ~ (|n — l\/2){Rh/ Rc)'^ , 
whose magnitude depends on halo mass, becoming smaller 
for halos of smaller mass. This is qualitatively in agreement 
with the results of A*'-body simulations discussed by Manera 
et al. (2010). 
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APPENDIX A: HALO-MASS CROSS 
CORRELATION WITH FOURIER-SPACE 
BIASING 

When the halo field in Fourier-space is linearly biased as 
in equation (|20p with 5pk — >■ 5h and _Rpk — >■ Rh, and with 
a constant bias b, then the bias determined from the ratio 
of the halo-mass cross correlation ghm and the mass auto- 
correlation function gmm, will in general be smaller than 
b. To see this, consider a toy model for the CDM power 
spectrum given by Ps{k) oc fce"'' ^c/^ Further assuming 
the halo field in equation (|20p to be Gaussian filtered with 



